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5y ABSTRACT 

We compute the divergent contributions to the one-loop action of the U(l) Abelian Higgs 
model. The calculation allows for a Friedmann-Lemaitre-Robertson- Walker space-time and a 
time-dependent expectation value for the scalar field. Treating the time-dependent masses as 
two-point interactions, we use the in-in formalism to compute the first, second and third order 
graphs that contribute quadratic and logarithmic divergences to the effective scalar action. 
Working in gauge we show that the result is gauge invariant upon using the equations of 
motion. 
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1 Introduction 



There is only one scalar field in the standard model, but it plays a crucial role. Scalars 
are more abundant in most extensions of the standard model, such as the multiple Higgs 
fields in grand unified theories, the superpartners in supersymmetric models, and the moduli 
fields in extra-dimensional set-ups. To study the physics of the early universe and to test 
physics beyond the standard model using cosmological data, it is important to have a precise 
understanding of the scalar field's dynamics. This requires to go beyond a classical treatment 
and include the dominant quantum effects. 

The one-loop effective action for a scalar field in an expanding universe has been known 
for a long time [U O [3l H] • It describes the backreaction of the quantum fluctuations of the 
scalar field on the (time-dependent) background field, which can be calculated systematically 
in a loop expansion. If the scalar is coupled to other scalars or to fermions via e.g. a Yukawa 
interaction, additional scalar and fermion loops contribute [2j[5j[6]. In this paper we extend 
these results by including a coupling to a gauge field. That is, we calculate the effective action 
for a Higgs-like field, which is charged under a gauge symmetry, in a Friedmann-Lemaitre- 
Robertson- Walker (FLRW) universe. For simplicity we focus on an Abelian symmetry, but the 
results are easily generalized to more general gauge groups. This is the generalization of the 
Coleman- Weinberg potential [7] to time-dependent background fields in a curved space-time. 

Our results are of direct importance for inflation models in which the SM Higgs or a GUT 
Higgs is the inflaton (or a waterfall) field (e.g. 019]). The Coleman- Weinberg potential gives 
the dominant quantum correction during inflation. The time-dependent corrections may 
become important at the end of inflation, and during the subsequent period of reheating. 
Another application is the description of flat directions of the MSSM and its extensions |10| . 
which are lifted by the one-loop quantum corrections. This may affect inflation models or 
Affleck-Dine baryogenesis models [11[ [T2] using flat directions. 

The effective action for an Abelian gauge theory in de Sitter space-time has been calculated 
by P3J HU H5] using the Landau gauge. More recently the calculation was done in the Rg 
gauge, showing gauge invariance of the effective action [16J. To obtain this result an adiabatic 
approximation was made which fails in the £ — > limit. We extend these results to a generic 
FLRW space-time and allow for the possibility of time-dependence of the background field, 
which in a cosmological set-up can be displaced from its potential minimum. The calculation 
is done in the R^ gauge using a perturbative approach. We calculate the quadratic and 
logarithmic divergent terms in the ultraviolet (UV) limit, which come from a finite number 
of diagrams (and which do not depend on the perturbation being small throughout). The 
resulting effective action is gauge invariant only on-shell, that is after using the classical 
equations of motion, in agreement with the Nielsen identities |17[ 118]. 

Our results agree with the expressions in the literature in the appropriate limit. In the 
limit of a static background field and a constant Hubble parameter our results agree with 
|16] , In the Minkowski limit we retrieve the effective action calculated in our previous work 
|19] , and also the effective equations of motion found earlier in [201 [2T1 [22] [23] . Finally, taking 
both a static background field and a static background we get the familiar Coleman- Weinberg 
potential [7]. 

To properly take into account the real-time evolution of the system, and to assure the 
results are manifestly real, we use the in-in formalism (also known as closed-time-path (CTP) 
or Schwinger-Keldysh formalism) to calculate the effective action [231 1251 1261 [271 1281 1291 [30] . 
We first derive the one-loop corrected equation of motion for the Higgs field, using the tadpole 
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method [31J. All one-loop Feynman-diagrams with one external Higgs leg contribute. For 
technical reasons it is easier to work in conformal time, as the resulting action has a form 
similar to the Minkowski action, and all the machinery developed for this [3"2] can be used. 
We split all two-point interactions into time-independent and time-dependent parts, treating 
the former as masses and the latter as interaction terms; in the loop-expansion the result will 
not depend on this split [TJ. The equations of motion can be formally integrated to obtain 
the effective action up to field-independent terms [33]. Finally, we can rewrite the results in 
coordinate time. 

The effective action is independent of the specific initial conditions chosen. We will argue 
that this is always the case, for arbitrary initial conditions, provided the initial vacuum is 
chosen to be that of the free theory. The different vacua can be related via a Bogoliubov 
transformation |34j . 

As already mentioned, we only calculate the UV divergent terms, as these will generically 
give the dominant contribution. Using a renormalization prescription, these terms (together 
with the wavefunction renormalization of the gauge field) suffice to derive the renormalization 
group equations (RGE) and find the RG improved action. We neglect the backreaction on 
space-time, and treat the FLRW scale factor as classical background field. Finally, we note 
that to apply the results to models of Higgs inflation, a non-minimal coupling to gravity has 
to be considered. All this is left however for future work. 

In the next subsection we give a self-contained summary of the results. Following this 
we go through the calculation, starting in Sec. [2] describing the model, the in-in formalism, 
and giving the vertices and propagators needed to compute Feynman diagrams. In Sec. Owe 
calculate the relevant graphs at first, second and third order which contribute to the one-loop 
equation of motion. These graphs are used in Sec. U]to compute the effective action for the 
charged scalar. Here we also present results when additional scalar and fermions run in the 
loop. Our choice of initial conditions, and their generalization, is discussed in Sec. [5j We 
conclude in Sec. EJ and provide some further details of the calculation in a pair of appendices. 



1.1 Summary of the results 



Here we shall outline the model and give the main result of our calculation. It is self-contained 
so that one need not get caught up in the details of the derivation to make use of the final 
answer. 

The ansatz for the space-time metric is 



ds z = dr - a z {t)dx 



(1) 



with a(t) the time-dependent scale factor of the FLRW metric. The action is that of an 
Abelian Higgs model, with a Higgs field charged under a U(l) gauge symmetry 



d 4 x, 



-g 



t 



(2) 



We expand the Higgs into a time-dependent background field (the zero-mode) <f> plus quantum 
fluctuations h and 9 



1 
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(4>(t) + h(t,x) + i9(t,x)). 



(3) 



We gauge fix using gauge (see details in next section). The effective potential is corrected 
by loops which are calculated using the in-in formalism. We sum over all the relevant loops, 
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and in the end go on-shell, which ensures the gauge invariance of the final result. We find, 
up to background field-independent terms, the UV divergent contributions at one loop to be 



r 



1— loop 



■1 



16tt 2 

where the "shifted scalar mass" is 



(V hh +V 0e +3m 2 A )A 2 - ( V% h + Vi e + 3m\ - W ee m\ ^ ' S ' "' ' 



V n , 



H -2H 2 , 



(4) 



(5) 



which is time-dependent. A subscript on V denotes a derivative with respect to that field. 
Further, A is the cutoff used in regulating the divergent momentum integrals (it cuts off 
3-momentum \k\ < A), and the arbitrary mass fh is put in to ensure the argument of the 
logarithm is dimensionless (recall that we are only interested in divergent contributions to 
the effective action). H = a/a is the Hubble constant, with a dot denoting a derivative with 
respect to time t. The (time-dependent) mass of the gauge field is m\ = g 2 (f> 2 ■ 

Our result agrees with those found in the literature. For the Minkowski case {H = H = 0, 
and thus V aa = V aa ) it matches our previous result |19j . In the de Sitter limit H = 0, and for 
a time- independent Higgs field (Vgg = by Goldstone's theorem), it agrees with Garbrecht 
[16] . Finally, taking both the Minkowski limit and a static background field, we retrieve the 
classic Coleman- Weinberg potential [7]. 

If the Higgs field couples to additional scalars Xa and/or fermion fields V^Q we get an 
additional contribution 



£pl-loop 



I 



16vH 



d 3 xdt- v /— g 



V ( V A 2 - V 2 HA/rh) 



Veem 2 ^ 



2 \ ln(A/m) 



(6) 



where the shifted scalar mass is given by (J5]). Here the sum is over all bosonic and fermion 
real degrees of freedom, where a Weyl (Dirac) fermion counts as 2 (4) degrees of freedom. 



2 Action and formalism 

The starting point for our calculation is the action of a U(l) Abelian-Higgs model in an 
FLRW background. The background space-time is fixed, in the sense that the backreaction 
of the charged scalar is assumed negligible. We work in gauge, and therefore include in the 
action a gauge fixing and a Faddeev-Popov term. We work with a conformal metric, where 
most expressions take a form reminiscent of the Minkowski calculation. Since the background 
space-time, as well as the scalar vacuum expectation value (its classical value), are taken to 
be time-dependent, the "masses" (really two-point interactions) of the particles are also time- 
dependent. We deal with this by splitting these two-point interactions into a time-independent 
part, which we call the mass and which determines the propagator, and a time-dependent 
part, which is treated as a proper two-point interaction in Feynman diagrams. The action, 

4 We assume here a Yukawa interaction oc <f>. The expression for a more general mass term is given by 
((57} . We work in a basis where both the mass and kinetic terms are diagonal; this can be easily generalized. 
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propagators and interaction vertices are denned in the following subsections, along with the 
in-in formalism for computing expectation values. We use this machinery in the next section 
to calculate the one-loop equation of motion. 



2.1 Notation 

We use a metric with signature (+,—,—,—), indexed by lower Greek letters fi,u..., and 
lower Latin letters for just the 3-space. The Greek letters a, /3, . . ., and I, J, index the set of 
quantum fields. For propagators, covariant derivatives, mode function normalization etc. we 
use the same conventions as Peskin and Schroeder [35J. Masses m and frequencies oj are split 
into time-independent and time-dependent parts, with the notation m 2 {t) = fh 2 + 5m 2 (t), 
with 5m 2 (0) = 0. A hat above a mass scale denotes the corresponding quantity rescaled by 
the scale factor: m = am. Derivatives with respect to conformal time r are denoted by a 
prime, and derivatives with respect to coordinate or physical time t by a dot. In momentum 
integrals dfc = dfc/(2-7r). 



2.2 The action in an FLRW background 

The FLRW metric in physical and conformal coordinates is, respectively, 
The non-zero connections are 

r i0 = r oi = = T° = %. (8) 

Here we defined % = a' /a, analogous to the usual definition in coordinate time H = a/a. We 
can decompose the charged scalar field into a real and imaginary part, 

$(^) = ^(^(r) +h(r,x)+ie(r,x)), (9) 

with </>(r) the time-dependent classical background field. The action is a sum of the kinetic 
and potential terms, the gauge fixing term and the Faddeev-Popov term: 

Stat = J d 4 x^(C + C GF + £ FP ), (10) 

with 



C = -\g^ a g^F^F a p + g^D^(D u $)t - V($), (11) 
£gf = -^-G 2 , G = g*'VrA v -Zg(<t> + h)e, (12) 

£fp = W-faV- ( 13 ) 

Note that S/ fl g^ 1 ' = (because of metric compatibility), and thus g^ v V ^A u = V^g' MU A u and 
there is no ambiguity. SG/5a is the operator obtained by computing the variation of G 
under a U(l) gauge transformation with infinitesimal parameter a. A bar over r] denotes the 
conjugate, not to be confused with a time-independent quantity. 
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Our aim is to compute the quantum corrected equation of motion for the background 
field 4>. It can be derived by demanding that all tadpole diagrams with one external h-leg 
vanish. To calculate these diagrams we need to derive both the propagators, which follow 
from the free and time-independent part of the action, and the interaction vertices, which 
include time-dependent two-point interactions, as well as three-point interactions!! The key 
to our approach, following [32], is the choice to separate out from the two-point terms the 
constant, diagonal pieces. This allows us to easily solve for the scalar and gauge propagators, 
since the masses are then really masses: they are field-diagonal and time independent. The 
left-over pieces are treated as two-point interaction vertices, and put into graphs following 
the usual Feynman rules. 



To make explicit all factors of the scale factor we now write g 



a rj' 



2 rf v with rf v the 



Minkowski metric, which is the metric in the comoving frame with coordinates (r,x). In the 
expressions below, all indices are raised and lowered using the Minkowski metric. We denote 
all mass scales in comoving coordinates with a hat. In particular we define: 



V = a 4 V(<f>), 



(14) 



with <j) a = {</>, h, 6, rj} the scalars in the theory. The hatted fields are canonically normalized 
in the comoving frame. Since the gauge field kinetic terms are conformally invariant, there 
is no rescaling of the gauge field. These comoving fields feel a potential V. All the comoving 
quantities map directly to the equivalent set-up in Minkowski, and we can use the usual 
Minkowski machinery to calculate Feynman diagrams. 

The action (110p is expanded in quantum fluctuations around the background. Here we 
state the results at each order; for details see Appendix [Aj The one-point vertex is ()104|) 



d 4 x 



-X h h 



where 



X h = (d 2 T -{W + U 2 ))^ + V^ 
The quadratic action is (|105j) 



+ 3H<p + 



(15) 



(16) 



5(2) 




<9 2 + m^rT + (1-7 ) QilQV 



A„-A (m 2 ) i0 A. 



m 2 Ae A t 



l - &(d 2 + ™lJ&-v(d 2 + m 2 )f)\. 



(17) 



<Pi={h,6} 

The explicit form of the two-point interactions, with rh 2 = a 2 m 2 , are: 



™00 



g 2 4> 2 + - {%' - 2V 2 ) 8^o = a 2 
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rh 2 =V hh -(W + U 2 ) = a 2 



V hh -{H + 2H 2 
mj = V ee + tg 2 tf - {%' + V 2 ) = a 2 



+ ^g 2 (/) 2 -{H + 2H 2 ) 



m 2 = ig 1 ^ 1 - {%' + n 2 ) = a 2 



2 a.2 



(H + 2H' 



(18) 



5 There are also four-point interactions but these do not contribute to the one-loop tadpoles. 
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where we used T~L 2 = a 2 H 2 and H' = a 2 (H + H 2 ). The off-diagonal two-point terms are: 

m 2 M = 2g{d T - H)j> = a 2 \2g$\ , (m 2 ) 0i = (m 2 ) i0 = %H& . (19) 

L J £ 

The mixing between the spatial Ai and temporal Aq gauge field contains a first derivative, 
which survives in all gauges except for unitary gauge £ — > oo. However, we will not work in 
unitary gauge as this is known to give false results if not done carefully (see for example |36j). 
and we are thus forced to deal with this extra complication. 

As stated above, we choose to split the two-point interactions into a time-independent 
and time-dependent part: 

rh 2 ( T ) = rh 2 + 5m 2 ( T ). (20) 

The first term contributes to the free Lagrangian from which the propagator is constructed, 
whereas the time-dependent term is treated as a two-point interaction. The loop expansion is 
independent of the split of the two-point terms into a free and interacting part [7] . The split 
is defined by requiring the interaction to vanish at the initial time, which we choose without 
loss of generality to be at to = 0: 

5m 2 (0) = 0. (21) 

For the diagonal two-point interactions, (|2ip serves only to define m 2 unambiguously, and 
does not constrain anything physically meaningful, like, for example, the initial conditions 
of the background scalar and scale factor. However, due to the presence of off-diagonal 
two-point interactions, and the fact that we want to retain a Minkowski-like propagator 
structure to simplify the calculation, we will need to make some choice. In particular, we 
choose initial conditions such that the initial off-diagonal two-point interactions and Lorentz 
violating contribution to the gauge boson mass (the term proportional to <5„o) vanish at the 
initial time: m^ ff _ 

diag(0) — This implies the initial conditions 

5(f)(0) = 5(f>'(o) =n(o) = n'(o) = o, (22) 

where we wrote (p(t) = <f> + 5(p(t). As we will argue in Sec. [5j the results we obtain do not 
depend on the specific initial conditions chosen, provided we define the initial vacuum as that 
of the free theory. It is therefore no real limitation that the boundary conditions above are 
not the most physically motivated ones in actual cosmological settings. 

The two-point interactions for the scalars, m^, m 2 , and m 2 , are split as per (|20p . The 
diagonal gauge boson two-point interaction is split as a constant degenerate piece, a time- 
dependent degenerate piece, and an extra piece for ^4° denoted by m^: 

^fti) = ™A + $ m A + <5"io ( W, (23) 

with 

Sm 2 = hn' -2U 2 ). (24) 

The off-diagonal two-point interactions do not have a time- independent part, so 5m 2 Ae = rh 2 A Q 
and (5m 2 )° l = (m 2 )°\ For these terms, we shall use the notation with and without the 5 
interchangeably. 

The three-point interaction vertices are (I106p 



= / d 4 x 



(25) 
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where ip a runs over the quantum fields {h,8,A}, and 



^hhh 


= d f l l = 


V<f>hh, 




= = 


%ee + 2*ig 2 4> 


^hAA 


= d ^A = 




A/ir^jj 


= d ^ = 




XhAO 


= M-9t 


-U). 



Note that XhA9 contains a derivative of conformal time, which, by using partial integration of 
(|lU6p . we let act on the factor Aq6 instead of h. This allows us to factor out a common h in 
and thereby compute the tadpole diagrams^ In ([26]) we express the diagonal couplings 
<^haa as derivatives, with respect to (p, of the corresponding two-point interaction defined by 
(fT8|) , In computing these derivatives, we have assumed a fixed FLRW background so that H 
is independent of <j) and d-/H' = diH 2 = 0. 



2.3 In-in formalism and propagators 

Since we are interested in expectation values of the background field, and their evolution 
with time, rather than scattering amplitudes, we will use the in-in formalism (also known 
as closed-time-path or Schwinger-Keldysh formalism) [211 [251 [261 [271 [281 [291 [30] . In contrast 
with the usual in-out formalism, this approach gives results that are manifestly real. 

Expectation values are computed using an action S = S[(pf] — S[(/)~], with boundary 
condition <fif(t) = 4>~(t). That is, we double the fields, and take the action for the plus- 
fields as for the minus-fields, given by the equations in the previous subsection. All fields, 
propagators and vertices are labeled by ± superscripts. In a Feynman diagram the propagator 
x — x') connects between a X ± (x) and a A ± (x') vertex. Since the action of the minus- 
fields is defined with an overall minus sign we have 

K /3 ]" = -K /3 ] + , K aP = -K aP - (27) 

We construct the propagators from the free, time-independent part of the quadratic action 
(|17p . The corresponding quadratic Lagrangian can be written in the form 

= -(1/2) Y, 0+(^)i^V)4(^), (28) 

a,f3 

with the sum over all fields (f) a = {h,9,r], A^ 1 }. For example, the scalar fields have K a ^ = 
(d 2 + rh 2 )S a P . The propagators are then defined as the solutions of 



K^(x^ \(D++{x»-y») D+-{x»-y») 

-K a P(x») ) I Dp+(x»-y») D^(xv-if) 



-iSZ6(x* - y*)I 2 . (29) 



This defines D ++ as the usual Feynman propagator, D as the anti-Feynman propagator, 
and D h and as Wightman functions. Our initial conditions fj21|22f) are such that the 



6 Equation (|106[) also contains a term —2gA'9dih. Since the final expression of each tadpole graph is 
independent of the spatial coordinates, this three-point interaction does not contribute to the overall result. 
We have checked this by explicit computation. 
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FLRW propagators in conformal coordinates are analogous to the usual Minkowski expres- 
sions. 

It turns out convenient to rewrite all the propagators in terms of Wightman functions. 
We introduce the shorthand for the Wightman function 

Di^ab = DJ + {x a - x b ), D^ ab = D^J~(x a - x b ), (30) 

for the propagator of a type-/ scalar, and the gauge boson propagator, respectively. Using 
this the propagators are (suppressing the I or the Lorentz indices) 

D ++ {x a - x b ) = D ab e ab + D ba e ba , 
d (x a - x b ) = D ab e ba + D ba e abl 

D~ + {x a - x b ) = D + ~(x b - x a ) = D ab , (31) 

with Q ab = B(r a — r b ) the usual step function. We Fourier transform the propagator with 
respect to comoving three-momentum 

D(x a -x b ) = j d?kD(r a - r b , k)e^~^\ (32) 

with (Ik = dk/(2ir). The time dependence (T a — T b ) in the Fourier propagator will from now on 
be suppressed. For a scalar field the solution for the Fourier transformed Wightman function 
is 

2uji(k) 

where Cdj(k) = fh 2 + k and fh 2 is the (conformal) mass-squared of the appropriate scalar. 
The Fourier transformed Wightman function for the gauge boson propagator is 

D„{k) = ~ Uu - k -^f] D A {k) - ^D^k), (34) 
\ m A J m £ 

with Da and D^ scalar propagators with mass-squared m 2 A and rn| = ^fh 2 A respectively. In 
the £ = 1 gauge the gauge boson propagator is diagonal, D^ = —t]^ u Da- 

The mixed two-point interaction (Srh 2 ) 01 contains a spatial derivative; it acts on the gauge 
boson propagator by 

(6m 2 ) 0i (T a )D^ ab (k) = -itt-HiTjDi^k), 

(5m 2 ) 0i (T a )D^ ab (k) = -(5m 2 ) 0l (n)D^ ab (k), (35) 

where we note that (Srh 2 ) 10 and D^ Ujab are diagonal in the Lorentz indices. The boundary 
conditions (|21|l22p imply that (m 2 ) 0i = fh\g = rfiQ = 0; the off-diagonal two-point interactions 
and the Lorentz violating Aq mass contribution only enter the interaction Lagrangian. 

3 One-loop equation of motion 

The equation of motion for the background Higgs field <j)(t) follows from the vanishing of the 
tadpole. In terms of diagrams, these are all one-particle irreducible tadpole graphs with one 
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external h + leg. In this section we compute these diagrams, and thus the quantum corrected 
equation of motion, at the one-loop level. We are here concerned only with the UV divergent 
contributions to the graphs, and thus need to consider only the diagrams with up to three 
vertices. The calculation is done in the conformal frame, in terms of hatted fields and mass 
scales, conformal time and momenta. For notational convenience, in this section we drop the 
hat on all quantities; it shall be reinstated at the end when we give the results. 



The calculation is analogous to the one for Minkowski [19], but with two-point interactions 

(|18p that now depend on the FLRW scale factor. This is straightforward to incorporate for the 
diagrams with a scalar running in the loop. There are however some new technical difficulties 
that come in with the gauge boson loops: 

1. The mass of the temporal gauge boson m\ gets FLRW corrections but the mass of 
the spatial components mf does not. This is possible because Lorentz symmetry is 
broken by the time-dependent background. Consequently the diagrams with Aq and 
contribute differently. 

2. The off-diagonal gauge boson two-point interaction (<5m 2 ) * is non-zero. This results in 
new diagrams with both two and three two-point insertions. 

3. The formalism is set up in such a way that the two-point interactions vanish at the 
initial time (|2ip . This avoids divergences that depend on the initial conditions. We will 
argue in Sec. [5] that this is always an allowed choice, for arbitrary initial conditions, 
provided the initial vacuum is chosen accordingly. 

The one-loop equation of motion can be extracted from the series of tadpole diagrams 
with one external h + leg, see Fig. [TJ This gives 



where Ai labels the i order contribution to all tadpole diagrams with i vertices. We labeled 
the classical contribution A c \, which comes from a tree-level diagram. Indices {oc,/3, . . .} are 
compound indices denoting field-type as well as possible Lorentz indices for the gauge field. 
The sum is over these compound indices (all fields and their Lorentz indices) as well as all 
possibilities for ±. The S a p... are appropriate symmetry factors, derived in Appendix|Bj 5m\g 
and Xh a p are the two- and three-point interaction vertices, respectively, as defined in Sec. 12.21 
The tree- level tadpole diagram contributes (|16p . and we recover the classical equations of 
motion (remember we dropped the hat for conformal coordinates and scales): 



Given the equation of motion, we want to find the corresponding action. We do this by simply 
writing down an action which, upon using the Euler-Lagrange equations, yields the equation 
of motion (I37D . This action is then transformed from the comoving to the physical frame, 




= A d + Ai + A 2 + A 3 + finite 




+ finite . 



(36) 




(37) 
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Figure 1: Tree- level tadpole giving the classical equation of motion and the first, second and 
third order diagrams respectively. The summation is over all fields, for the gauge bosons also 
over Lorentz indices, and over ± at the two-point vertices. 



thereby obtaining the effective action. This will be done in Sec. [H The idea is to apply this 
procedure to the quantum corrected equations of motion in order to determine the quantum 
corrected effective action. 

We divide the calculation based on the order of the contributing graphs, which is the 
number of vertices in the loop of the tadpole. As discussed above, we must work to third 
order. Independent of this, we can distinguish three classes of diagrams depending on how 
they contribute to the answer. First there is the contribution that is fully analogous to the 
Minkowski calculation A Mink = Af mk + Af' mk , the only difference is that the mass term of 
the scalars now depends on the scale factor. Second is ^4 mass = ^4^ iass , which arises from the 
extra Feynman diagrams due to the FLRW mass correction of the temporal gauge field 5m^; 
see (|23p . And finally yl mix = A™ lx + A™ x gives the diagrams with one and two off-diagonal 
vertices (Sm 2 ) " 1 connecting the temporal and spatial gauge fields, also absent in Minkowski. 

3.1 First order contribution A\ 

The calculation of the first order diagrams proceeds analogously to the equivalent calculation 
in Minkowski, which can be found in |19| . At first order, four diagrams contribute, with 
ip a = {h,6,r], A^"} running in the loop. The result only depends on the time-independent part 
of the two-point interaction, as there is no vertex insertion. For each diagram the result has 
the same structure, given by (|113|) 

(Af nk ) a = ~^m^++(0). (38) 

Just as in Minkowski the gauge loop can be expressed in terms of scalar propagators (|34p via 

- rT £+ + (0) = 3D++(0) + £Z>++(0). (39) 



II 



The sum of all first-order diagrams is 



A 



k 2 dk 



1 mi 



+ ... 



16-7T 



r-2 ^ fi^S, 



1 



A" 5 m z a \n(A/m) 2 + finite 



(40) 



with k < A a cutoff on comoving momentum. The sum is over a = {h,6,r],A,S,}, and 
$a = {1, 1, — 2, 3, 1} counting the real degrees of freedom (with a minus sign for the anti- 
commuting ghost). Further m 2 A = g 2 4> 2 and m| = £,m 2 A . Note that the factor d a m 2 a is 
time-dependent, and evaluated at r; hence ^4^ mk is a function of r. The finite terms that we 
have neglected remain finite as A — > oo. 



3.2 Second order contribution A 2 

At second order the loop diagrams with one two-point insertion contribute. We split them into 
three parts. The first part, A% contains all scalar loops, and the gauge boson loop where 
only the diagonal part of (|23|) is inserted. In addition there is a mixed 9 A -loop (|115p . This 
part is analogous to the equivalent Minkowski calculation. The second part yl 2 nass contains 
the gauge boson loop with a Stuq, and the third part A™ 1X , with a (<5m 2 )°\ These last two 
diagrams are both absent in Minkowski. 



3.2.1 Af ink 



The scalar Higgs loop with one two-point insertion gives (|114|) 



a Mink 
A 2,h 



--d^m^Ta) I d x b 5m h (T b ) D+ + (x a -x b )D+ + (x b -x a ) 
- D^~(x a -x b )D^ + (x b -x a ) 



d^m h {r a ) / d x b 5m h (T b )Q ab 



Dh,ab ~ D 2 ba 



(41) 



where we expressed the result in terms of Wightman functions, using the notation introduced 
in Sec. 12.31 Fourier transforming as in (|32|) and performing the d 3 x integral gives a 5 s (k + p). 
And thus: 



a Mink 
A 2,h 



d$m%(T a ) / dT b 5m h (T b ) / 3 k 



d<pm h {T a ) I dT b 5m h (r b ) 



d 3 k 



Dh,ab{k)D h>ab (p) - D hM (k)D hM (p) 
sin [20 h (T a - T b )] , 



k=—p 

(42) 



(2^) 2 

where we used the explicit form of the propagator (i33j) . Now use integration by parts to 
extract the UV divergent piece (we give here the general formula, which can also be applied 
to the gauge boson loop discussed below) 



Jo dnf{Th) J (2u>r)(2u>j) Sin[( ^ + ° j){Ta " U)] = f(Ta) J ^ 



1 



(2w 7 )(2wj)(w7 + wj) 



(43) 
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where we assumed /(O) = and the ellipses denote higher order terms Oik 4 ). Putting it all 
back together we find for the scalar loop 



Aft = VkW) / = -d^mi{r)5mi(r) 



d 3 fc „ n , N „ n , , f d 3 fc 



(2^) 3 v n " ' " v ' J 8P 
= -d+ml {r)6m 2 h (r)— ^ m(A/m) 2 + finite, (44) 

where we rewrote the time variable r a = r. To get to the second expression we expanded in 
large momentum. 

The Goldstone boson and ghost loops give a similar contribution, although for the ghost 
with an overall factor (—2) to take into account that these are two real anti-commuting degrees 
of freedom. The calculation of the gauge boson loop follows the same steps, except that now 
care has to be taken of the Lorentz structure. We find 

A Mink = _^ m 2 J d A Xb6m 2 ATb)rj ^ vP , D ++(x a -X b )D++(x b -X a ) 

= -d^m 2 A \ dT b 5m 2 A (n) / 3 k sin [(a)/ + uj){r a - r h )\ 

Jo J \*ui){2u)j) 



= -d^m\5m 2 A {T) ^- \ 1 ln(A/m) 2 + finite. (45) 
Again, we rewrote r a = r. The relevant propagator combination, which defines C/j, is 
rr<rD lv (£)D m @)\ s __ / =Cij(k)D l D J (46) 

= f 3 + ^P] D A (i? + (' ( 1 + DM - 2 ^% + ^ 2 D A (k)B <( i), 

with I, J = A,£, and D Jia6 = (aS/)-^-* ^ 7 ""- 7 *). 
Finally, the mixed #j4°-loop gives f j 1 1 5 [) 



45$ = -*A WA (r a ) / d 4 x b <5m Ae (T b ) 



r>++ n++ _ D +- D -+ 

^ 00, ab^efia ^ 00 ,ab^ 6 ,6a 



= -2A WA (r a ) / dr b 5m 2 Ae {T b ) / 3 fcV . w r sin [(a)/ + a) e )(r a - t 6 )] 
Jo J (2wr)(2u; j 

= 2A MA (r)*m^(r)^^ ln(A/m) 2 + finite, (47) 

with A/ 1 6»a(t) = 2g(—d T — %{t)) the appropriate three-point vertexQ In the second line we 
identified 

(48) 

In the final step we performed integration by parts (|43p . and took the large k limit. 



7 Note that J &T\h6A(T)8m 2 Ag {T) = J drd^m\ g (T)Sm, Ag (T) after integration by parts, which we will use in 
the expression for the effective action in Sec.[J]to rewrite {A^ lnk )Ae in the same form as the other contributions. 
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Adding everything together gives 

= "3^2 E Sad^mlSml In A 2 + ^±^-\ hAe 8m 2 M ln(A/m) 2 + finite. (49) 



3.2.2 Af &ss and Af lx 

The last two diagrams to contribute at second order are those with an m§ and an (m 2 )° l 
insertion, the Lorentz violating mass and off-diagonal gauge boson interaction respectively. 
Neither diagram is present in Minkowski. The loop with m 2 , gives (|114p 

^ asS = -\d,m\(r a ) I d 4 x 6 5m 2 (r b )^ [d^U^M ~ rfo^M 



-d$m 2 A (r a ) J dn5ml(T b ) J d 3 /c^ 



Cu 



sin[(w/ +cjj)(t - n)] 



= -a ( ^m 2 4 (r)5m 2 (r)^±0 l n ( A /m) 2 + finite. (50) 
Here we defined the relevant propagator combination by 

VV"»m>(£)»o,®|j,_ ? = ^C7jZM£)A;M (51) 

m A m| m A m^ 

with as before I,J = A,£, and Dj(-jfe) = £>/(£) = l/(2w i )e- io '( T -" 7 *). 

The off-diagonal interaction (<5m 2 ) * contains a spatial derivative, and brings down a factor 
of the momentum. The diagram is given by (|116[) 

r,++ r,++ _ D +- D -+ 

fi0,ab iu,ba fj,0,ab iu,ba 



- k=—p 



Af* = id^m 2 A {r a ) J d 4 x b {5m 2 t(r b )t v 

2 r T a r 3 

= -|^m A (r a ) y dT b H(T b ) J ft kp 1 ^" D abttl0 (k)D abiiu (p) + D ba4l0 (k)D batilJ (p) 

= -^dsm 2 A {T a ) / dr 6 ^(r 6 ) / d k 2cos [(wj + wj)(r a - r b )] 

? Jo J {2ui){2ujj) 

2 P> 2, ^ />, 2ft'(r a ) 
= -~d< j> m A (T a ) / d fc _ Cjj + finite 

WCtVI — z^ 2 

= 3^(r) " ^ ^ ln(A/m) 2 + finite. (52) 

As we now have a cosine instead of a sine in the expression on the third line above, we 
integrate by parts twice to isolate the leading term in the UV limit. This is why the result is 
proportional to W . The relevant propagator contribution is defined by 



k 2 cu A (2k 2 + fh 2 A ) 2 k-u&k'+ml) k 2 {u A +00^ +oo A OJ i ) 
Tl L>A ^4 D £ H ^4 L>AL>i- 

m A rriA m a 
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3.3 Third order contribution A* 



The third order diagrams with two two-point insertions are UV finite, which can be easily 
checked by power counting. The only exception to this is the diagram with two off-diagonal 
(m 2 )° l insertions, because each insertion contains a spatial derivative, and thus brings down 
a power of momentum. We thus consider the third order diagram with two mixed-interaction 
insertions (|117p 



p. a 



where the sum a is over the four possibilities for the Lorentz indices 
(p, a, k, t) = (i, 0, j, 0), (0, i, 0, j), (0, 0), (i, 0, 0, j) 



(53) 



(54) 



and also over the four possibilities for plus and minus fields. The x a vertex is a +-vertex. 
There are then 4 possibilities for the vertices: 



(a,b,c) = (+ + +), (+-+), (+ + 



+ 



(55) 



The propagator between a (±)-vertex and a (rt)-vertex is -£> ±=l= , and we write them out in 
terms of Fourier transformed Wightman functions using the notation of section (|2.3p . Taking 
the action of the spatial derivatives in (m 2 ) * will then bring down powers of momentum, as per 
(|35p . Finally, the x b and x c integrals give delta functions encoding momentum conservation. 
At the end of the day we find 



o 



Af x = -d^m^Ta) I dr b j dT c U{T b )U{T c ) j d 3 ^^, 







(56) 



D ba {k)D bc (p)D ac (q) + c.c 



k=q=—p 



2Q bc (D ab {k)D bc {p)D ac {q) + c.c 



k=p=—q 



where the sum is now only over the Lorentz indices (|54j) , which we suppressed in the above 
formula. The sign s p = (1, 1,-1,-1) for the four possibilities ()54p . The relevant propagator 
combinations, putting Lorentz indices back in, are 



Smk'V D ba {k)D bc {p)D ac {q) 



k=q=—p 



k l Wrf v D^(k)D 0j (p)D 0u (q) + D^(k)D i0 (p)D JU (q) 

D^DijipjDo^q) - D pi {k)D m (v)D JV {q) 
Y,CiJK(k)Dj(k)Dj(k)D K (k), 



k=q=—p 

(57) 



and 



Smk'k? D ab {k)D bc {p)DaM 



k=p=—q 



: k l Wrf v D^(k)D 0j (p)D 0u (q) + D^{k)D iQ {p)D JU {q) 

- k=p=—q 

: ' S ^D I j K (k)D I (k)D J (k)D K (k), (58) 
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with I,J,K = A,£, and C IJK ,D IJK ~ k 2 . Using D I;Clb (k) = (2cd I )'- 1 e~ iiil( - Ta ~ T ^ we have 



A 3 -72 



V&fa.) jT° dr fe jT° dT c H(n)U(T c ) J d 3 k 



{2u i ){2uj j ){2uj k ) 
Cijk cos ((ui K - w/)r a + (a)/ + wj)r b - (wj + w^)r c ) 



- 2Q bc D IJK cos ((w/ + w^)T a - (uji - ujj)T b - (uj + wk)t c ) 

Now use integration by parts with respect to r b and r c to write 

dT b U(n) / dT c U{T c ) cos [(a)/ + a)j)r fc - (caj + a)x)r c + (% - w/)r a ] 
Jo 

where we used the initial conditions (|22p . And similarly 



/■Ta /"T a 

/ dT b H(T b ) dT c n{T c )@ bc COs[(idl +UK)T a - (tOl - Oj)T b - (UJ +U) K )t c 

JO Jo 

fr a Win) 2 

' dT b — — _ sin [(a)/ + w^)r a - (a)/ + u)k)u] 

(Wj + OJK) 

■H(r a ) 2 



(ujj + uj k )(uj i +L0 K )' 



The end result is 



/imix 

^3a - 72 



a^mi(r)^ 2 (r) | d 3 k^ 



1 



(2wi)(2wj)(2wj<-) 



0$m A {rjH. (r) ln(A/m) + finite 



64vr 2 £ 



64vr 2 £ 2 
ln(A/m) 2 + finite. 



3.4 Summary of graphs 



(59) 



(60) 



(61) 



(62) 



In the previous subsections we have computed all quadratically and logarithmically divergent 
contributions to the one-loop equation of motion. Here we collect and summarize the results, 
putting the hats back on the relevant variables to indicate that we are still in the conformal 
frame. 

The first order graphs are given by (|40p . The second order contributions are ([49|5QI52p . 
and are summarized in Fig. [2 At third order there is only one piece, given by (162j) . We now 
collect these terms into the three groups A Mmk , ^4 mass an d A imx . 

is 



The first and second order combined A Mink = Af ink + Af ink 



{Mink 



167T 



A 2 --m 2 lnA 2 

2 a 



+ ^rXhAemle ln(A/rh) 2 . (63) 
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+ 




-\hAe5rh 2 M ^ 



Figure 2: The second order tadpole diagrams and their corresponding mathematical expres- 
sion (below each graph). These Feynman diagrams are in (conformal) coordinate space, with 
the left and right vertices at x a and Xb respectively. The argument of each of the propagators 
is (xb — x a ), and all time-dependent quantities (Xao, Thai anc ^ ^) are evaluated at r. 



As expected, this is independent of how the two-point interaction is split into a free and 
interacting term, since the first and second order pieces combined in the sum = m^-\-5m^. 
For Aq mass insertions we have the second order piece (|5U|) 

i-ss = 5lh g3+j£ ln{k/fh) 2_ (64) 

For the mixed piece we have contributions from second order (|52h and third order (162p , giving 
a total 



A- = d^ A ^ ^ -I ^ ) 64^ ln ( A /-) • ( 65 ) 

All factors in ([551 ES ES]) that are time-dependent — being the m 2 's, and H — are 

understood to be evaluated at r. 
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4 Effective action 



The previous section found the one-loop equation of motion. The corresponding loop-corrected 
effective action is the one which, upon applying the Euler-Lagrange equations for 0, yields 
the loop-corrected equation of motion found in Sec. [3J The effective action is defined this 
way up to an arbitrary field- independent constant, and an overall minus sign which is fixed 
to obtain the correct sign kinetic terms. Our working assumption is that the background is 
fixed, i.e. a(r) not a function of the background field 4>(t). 
The classical action is defined by 



r 



cl 



d 3 x dr C 



cl 



where A d is the equation of motion following from the Lagrangian density £ cl : 



A c 



5C C 



5C C 



(66) 



(67) 



From (1371) we have that A d = A + , so the classical action is 



d^xdr 



1 



V 



d 3 xdr v / -5 conf 



J d 3 xdt v / -5p hys 



a 

1 

'2^2- 

-^{d 2 + md t )^-v 



v 



(68) 



where we used that %' + T-L 2 = a" /a, and H = a/a. In the second line we went to unhatted 
quantities, and in the third we changed to physical time. The measure in conformal coordi- 
a 4 , and in physical coordinates yj 



5ph y s = or- Taking the Euler-Lagrange 
we get the familiar FLRW equation of 



nates is y/-g CO ni 

equation from the last line, using the definition ([6 
motion: <p + 3H<p + = 0. 

Applying the same procedure to the quantum corrected equation of motion gives us the 
quantum corrected effective action. The relevant terms at the level of the equations of motion 
are summarized in Sec. 13.41 and the one- loop correction to the effective action is defined as 



pl-loop = J d 3 xdr (£Mink + ^mass + ^mix + fafay 



(69) 



In determining C, we must be careful to use the same sign convention in the Euler-Lagrange 
equation as we did above in (f67|) . to ensure the corrections to the effective potential have the 
correct relative sign. 

All but one term in J 4 Mmk are polynomial, the exception being the XhAe term. For this 
term, the <fr dependent factors are 

\hAem 2 A = V + n'4> + U 2 $) ■ (70) 
This expression comes from a Lagrangian 



: m Ae 



i/2 



2 22 



(71) 
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For the rest of the terms in ^4 Mmk 5 which are polynomial in cf>, the corresponding action is 
found simply by integrating with respect to </>, and then negating (since the 5C/8(j) term in 
(|67p comes with a minus sign). All terms in A mass and A mix are also polynomial in 0, so can 
be similarly integrated. Thus, from ([631 [Ml [65]) . and using (|7T|) . we obtain 

OL / " 

' mass - 1 ,m\ Mk/fhf ( (3+ / 2) (2^ 2 - %'] 



64tt 2 a k ' ' \ i 
= ~K ln(AM) 2 (^-J?H' - ^±^ 2 ) , (72) 

with a = {h, 9, rj, A, £}, and 5 a = {1, 1, —2, 3, 1}. 

Now take out a factor y/—g, and write the hatted variables in terms of their unhatted 
counterparts to go back to the physical frame. Use that Ti 2 = a 2 H 2 and H' = a 2 (H + H 2 ), 
and transform to physical time. All terms are proportional to the fourth power of a mass, 
hence we factor out an a 4 . The result is 

T 1 - 100 ? = J d 3 xdt^g~ (£ Mink + £ mass + £ mix + finite) , (73) 

with 

£ Mmk = -T^E 5 « (< A2 ~ \< ln(A/m) 2 ) - ln(A/m) 2 , (74) 

W"4MAM) 2 (^ 

C °"" = - 64^"* ln(AM)2 (^f^<" + ^) - ^^H 2 ) ■ (76) 

Whereas A was a cutoff on conformal three-momentum (equivalent to comoving momentum), 
A is now a cutoff on physical three-momentum. 

When we now plug in the FLRW-corrected two-point interactions (|18p . we find for the 
total one-loop effective action (up to field-independent terms) 



£ mass = _ 1 m 2 AHA/fh) 2 ( ^±^1 {H 2 - H) ) , (75) 



[V hh + V ee + 3mi) A 2 



1 '/,,, - // -111 1 ) 2 +(v ee -H -2H 2 ) 2 + 3m\ 
+ 2£V ee m A - (6 + 2i)g 2 ^ 2 + 6m 2 A (H + 2H 2 ^ ln ( A /™' 2 



) 



4 



(77) 



Recall that m 2 A = g 2 4> 2 . This result is still gauge variant, which was to be expected. Gauge 
invariance is only achieved on-shell. For a time-independent situation (4>(t) = const.) the 
Nielsen identities [IB] 

dV e s dd> dV c s , 
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show that the effective potential is only gauge invariant when the background field is in 
a minimum of the potential. Here (just like in [19]) we want to use the time-dependent 
version of this statement: the effective potential is only gauge invariant when the background 
field satisfies its equation of motion. Going on-shell enables us to rewrite in ()77p the term 
proportional to (ft 2 . This term originated from the mixed Goldstone-gauge boson "mass", the 
last term in (|74p . and can be transformed to: 



j4 / — 4 
d xy/-gm Ae 



d^x^Ag 2 ^ 



d^xyj— g4m A ] 



(79) 



In the third step we integrate by parts and go on-shell. The last step uses Goldstone's theorem 
(it exploits the fact that the potential is a function of [371 [Ml US])- 
On-shell the result (|77p takes the form 



-il— loop 



16vr 2 



d 3 xdt^/^g 



(V hh + Voo + 3m 2 A )A 2 



V hh -H-2H 2 ) + 



H-2H 2 )' + 3m 4 f 



, m \(v ee -H-2H 2 )y^M 



which is gauge invariant, as it should be. Introducing the notation 

Vaa = V aa — H — 2H 2 , 
we rewrite the final result (up to field-independent terms) 
-1 



(80) 



(81) 



pi— loop . 



167T 2 



d• i xdt^/^g 



{V hh +V ee +3m A )K 2 -[ V£ h + V e 2 e + 3m\ - 6V ee m 2 A 



ln(A/m) 



(82) 



4.1 Fermions and additional scalars 

It is straightforward to add fermions and additional scalars to the calculation. If these fields 
are coupled to the Higgs field, and thus have a ^-dependent mass term, they will contribute 
to the effective equation of motion for the background Higgs field <f>{t) and to the effective 
action. 

We assume the extra scalars are in a basis with canonical kinetic terms and have diagonal 
masses, and do not mix with h. Similarly, we assume the extra fermions have diagonal masses. 
It is easy to relax these assumptions and generalize the results. 

In terms of Feynman diagrams, there are extra tadpole graphs with the additional scalars 
and fermions running in the loop. The calculation for additional scalars is analogous to that 
of the Higgs fluctuations h already done, with a contribution at first and second order. The 
result is 



il— loop 
(scalar) 



1 



16vr 2 



d?xdtyf^g 



XX 



H - 2iT 



ln(A/m) 



(83) 
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where x is the additional real scalar and V(x, 4>) its potential. 

Just as for the bosons, the tadpole diagrams with a fermion loop can be mapped to the 
calculation for Minkowski space, except that the "mass" terms now depend on the FLRW 
scale factor. To discuss fermions in curved space-time, one has to use the vielbein formalism 
to transform to a local Lorentz frame, where Lorentz transformations and spin-^ particles are 
well defined. The vielbeins are defined via 



fj t n 



(84) 



with ef. 

H 1 



a<5" for a conformal FLRW metric 0. The gamma matrices are {7^, 7 1 '} = 2g< 



with 7" = e^7^ the usual Minkowski gamma-matrices. With this notation the fermionic 
action is El 



C f = / dW^M^V,, - m)</>, 



(85) 



with the covariant derivative V, 



+ Q^i, and f2 



(l/4)w a6M 7 a 7 6 , with Q = and 
0^ = (l/2)(a'/a)7*7° for the conformal FLRW metrical We rescale the fermion field ijj = a 3 ^ 2 ip 
and mass rh^ = am^. The Dirac equation then becomes 



{vfdp - rhip) %j> = 0, 



(86) 



which is of the usual Minkowski form. Hence the end result is the Minkowski (Tj result but 
with the replacement — > rh^ = am^ [6]: 



-il— loop 
(fermion) 



f J 



™W ~ 4 



™% + rn'^rh^) ln(A/m) 



16vr 2 



m\k 2 



fh^ + 3Hrh 



+ H + 2H 2 \ ln(A/m 



(87) 



The sum is over all fermionic degrees of freedom, which are two (helicity) states for a Weyl 
fermion and four states for a Dirac fermion. The first line is the Minkowski result with the 
replacement — > fh^. In the second line we went to physical coordinates by factoring out 
an overall a 4 factor, and rewriting the rh'^ in terms of derivatives with respect to physical time 
t. The first contribution to the logarithmic term incorporates the expansion of the universe. 
The second contribution to the logarithmic term is because the cj) field is rolling, and is also 
present in Minkowski space-time. 

Again we can simplify this result by going on-shell. For a fermion mass = X(p that 
is linear in the Higgs field — which is the case for Yukawa interactions and also for gaugino 
masses in supersymmetric theories — this gives 



il— loop 1 

(fermion) 167T 2 



E 



m\k 2 



- m^Vggj ln(A/m 



Here we have used, again, the background field equations and Goldstone's theorem. Vgg was 
defined in (l8il. 



If the fermions are charged under gauge groups, there will be an additional gauge connection. These extra 
terms do not affect the effective action for <f>, and for simplicity we leave them out. 
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5 Initial conditions 



Our interactions are time-dependent, and thus we needed to define the split between a time- 
independent mass and a time-dependent two-point interaction (|2ip 

m 2 aP {t) = mlp + Smlpit), 6m 2 a p(0) = 0. (89) 

We furthermore chose initial conditions for <p(t) and a(t) such that the off-diagonal and 
Lorentz violating two-point interactions vanished completely at the initial time (|22p . These 
choices ensured the simplicity of the propagators. They also ensured the vanishing of the 
t = boundary terms coming from integration by parts when evaluating the loop diagrams 
in Sec. El If these boundary terms did not vanish, they would yield extra contributions to the 
final result, contributions that depend on the initial conditions, and that diverge as t — >• 0. 

Our chosen initial conditions are peculiar, and are not the ones to be used in a realistic 
situation. The problem in straightforwardly generalizing our calculation to arbitrary initial 
conditions are the two-point interactions (m 2 ) *, m^, m^. To simplify the structure of the 
free action, and use the standard expressions for the propagator, we have treated them as 
interactions m 2 ^ = Srrv^g. To satisfy (j89j) then requires the initial conditions (|22l) . 

However, in principle there is nothing to stop us from also splitting these two-point in- 
teractions into a free and interacting part, as in (189(1 . Technically, this is complicated, as 
Lorentz symmetry is broken, and the gauge fields and Goldstone bosons all mix at the initial 
time. Nevertheless, in principle we can expand all fields in mode functions, where the mode 
functions satisfy the off-diagonal mode equations (diagonalizing the equations will result in 
a momentum-dependent diagonalization). Then (|89p is satisfied, all terms depending on the 
initial conditions vanish, and the results are the same as for our choice of initial conditions 
d22|). 

In slightly different words, we argue that the result is independent of the initial conditions 
as long as we choose the initial vacuum to be that of the free theory, which is defined by the 
split of the quadratic term into a time-independent mass and a time-dependent interaction 
term. That is, solve the mode equations derived from the free action with fh 2 ^^, and the cor- 
responding annihilation operators annihilate the vacuum. The different vacua, corresponding 
to different initial conditions, are then related by a Bogoliubov transformation. For the scalar 
field theory this was shown by Baacke et al. |34] (see also [39]). For U(l) model we require a 
more general Bogoliubov transformation, with momentum and polarization dependent coeffi- 
cients that mix the fields. In principle this should be straightforward, but we will not present 
any further details here. 

In practice, choosing the initial conditions ()22[) simplifies the calculation of the free field 
mode functions and propagators, and eliminates boundary terms, which is why we choose it. 
We have argued that a full treatment of initial conditions would yield the same result, at least 
for the divergent corrections to the equation of motion. 

6 Conclusion 

In this paper we have computed the one-loop divergent corrections to the effective action of 
a U(l) charged scalar, whose background vacuum expectation value is changing with time, 
and when the background space-time is of FLRW form. We used the in-in formalism and 
gauge, and our main aim was to demonstrate gauge invariance of the one-loop corrections. 
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The gauge invariance is only manifest upon using the equations of motion, i.e. on-shell, in 
accordance with the Nielsen identities. Our result is given by (|82p . and directly generalizes a 
previous work [115]. In comparison with the result in a Minkowski background, one can obtain 
the FLRW correction by shifting all scalar masses by 2H 2 + H. We showed that additional 
scalars running in the loop can be easily accommodated (|83|) . a result which also mimics the 
Minkowski case, but with a shifted mass. For fermions in the loop, the additional correction 
is (187p . Although our assumptions for the initial conditions of the background scalar and 
scale factor are unrealistic, we argued that, if handled in a more general way, the result would 
be the same. 

The computed correction can now be used to derive the renormalization group and find 
the RG improved action, a task we defer to a future publication. An additional task left for 
the future is to take into full account the backreaction of the scalar on space-time. Essentially, 
one must allow for spin-0 fluctuations of the metric, determine their mixing with the scalar, 
diagonalize to a new basis, and use this basis as the starting point of the calculation. A further 
generalization is to include a non-minimal coupling to gravity, so as to describe models of 
Higgs inflation. Finally, one could also generalize the decomposition of <1> ([9|) to allow for a 
time-dependent classical background in the imaginary direction. 
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A Action in detail 



Here we work out the explicit form of the action (jlOp to fourth order in quantum fluctuations. 
Using conformal coordinates the overall volume factor is y/—g = a 4 , and g^ u = a~ 2 r]^ u . Unless 
otherwise stated, all indices below are raised and lowered using the Minkowski metric. 

Start with the kinetic term for the gauge field. The connection cancels in the field strength 
Fa V = VuA„ - V v An = dnA u - d v Aa, and thus 



(90) 



As expected, the result is invariant under a conformal transformation of the metric. The 
kinetic terms and potential for the Higgs field are expanded as 



d 4 x-/=5(|£>$| 2 



V) 



J d 4 x| y Yl ((^) 2 + ^V) + 2gA»(-ed^ R + fad. 



<P=<PR,b 



d 4 x 




a 4 U 



+ gA» 



7; ( Vid 2 ~ ^ + a 2 V w )ip - g 2 A 2 (p 2 ) - ^aV^if - ^U ww <^ 4 



1 

3! 



-aOd, 




(91) 



23 



with 6r = (f> + h. The prime denotes derivative with respect to conformal time r. We rescaled 
the scalar fields tp a = {6, h, 9} as in (fT4"j) , 
The gauge fixing action is 



Sgf = ~ J d 4 *^ [OT V^) 2 - (g^V,A„)2£g6 R 9 + fg 2 ^ 



(92) 



The first term becomes 



-1 



-A^d u A v + A {%' - 2U 2 ) A - A 2Hd i A i 



(93) 



To get the second line we used the explicit form of the connections ([8]) to write 

Vrr£, = A ( V 00 T 00 + rj^Tl) = A U(l - 3) = -27iA , (94) 
and integration by parts 

4 J d 4 xA Q HdoA = -2 J d A xA n'A . (95) 

The second term in (|92p we can partially integrate using (with A* 1 = g^ v A v to indicate 
the index is raised with <T ) 



^x^g-(y^)B = - / d^y^i^S. 



(96) 



This follows from the fact that we have a covariant volume and a covariant derivative. Note 
also that V ' u.9u,v = 0, and it is therefore irrelevant whether the raised index is on A or on V. 
Thus the second term in (1921) can be written as 



d A xgA^ 



aOd,, 



+ a6 R d. 



a 



(97) 



The second term above will cancel with the last term in (|9ip . The complete gauge- fixing term 
is 



Sgf = J d 4 x 
-gA» 

Finally the Faddeev-Popov term is 



A^d v A v + A (W - 2V 2 ) A - A 2Ud i A i 



<>[d»- ^1) 4>r + 4>r (fy - ) 



S FP = J d 4 xa 4 f] [- V 2 + £g 



' - <Ph) I »7, 



which follows from 



5 a G = --3 2 a - -T^dPa + ig{6 2 - 6%)a = --V 2 a + £g(6 2 - 6 2 R )a, 
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9 



9 



(98) 



(99) 



(100) 
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where we used 5 a <f>R = —a9, 5 a 6 = a(f>u, and 5 a A^ = (—1/ g)d fl a. Use T pp = d p ^J —gj —g 
to write the first term in (1991) as 



d x^—gr]X7 r] 



d X\/—gr] 



1 



-(] 





1 d 4 xfi 


y - °- 


V 






a 





where in the last step we rescaled the anti-commuting scalars r] = arj. Hence 



Sfp 



I 



d 4 xr] 



d 2 - u - + ^l-e 2 ) 



fj, (101) 



(102) 



Putting it all together, we write the action as S = & with % denoting the number of 
quantum fields each term in contains. Then 



5 (o) 
5(2) 



<1 : Al^) 2 + ~< 
I a 



d?x{ -h[(d z - —)<p + a 6 V 4 



(103) 
(104) 



(d p dP + g 2 4> 2 )r]^ - (1 - hd^d" 

4 



+ - [A {%' - 2H 2 ) A - A 2Ud l A l ] 

- 9{d 2 -- + a 2 V ee + ^ 2 )6 - AgA°9 (d T - °L 
a \ a 



h(d 2 



+ a 2 V hh )h - 2fj 



S 0) 



d 2 -- + ig 
a 



212 



(105) 



S a 1 aV a i3 1 (pa<fl30'y 



2gA»§ (d„ - ^S°^j h + g 2 (A 2 - i§ 2 - 2#?)# j , 



(106) 



S {4) = j d 4 J -Saf hS Vaf h 60a<P{}<P'y<PS + \g 2 A 2 (h 2 + 2 ) - g 2 £fj(0 2 - h 2 )fj - ^faPh 2 \, 



with ip a = {h,9}, and S a p^(s) symmetry factors. For a quartic Higgs potential 



a 4 V = a 4 



^ 2-2,^/4 



1 a 2 22 i ^ 24 



(107) 



(108) 



B Tadpole method 

In this appendix we give the derivation of the one-loop quantum corrected equation of motion, 
starting from first principles to determine the symmetry factors of the diagrams. We use 
the in-in formalism (see Sec. Sand [231 E51 123 E3 129 ESI GS3) where all quantum fields 
are doubled, and labeled by ± superscripts. Expanding the Higgs field around the classical 
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background as in ([3]) , the the quantum corrected equations of motion follow from the vanishing 
of the tadpole [30] : 

{h+{r,x)) = 0, (109) 

where < • > denotes the vacuum expectation value. The vanishing of the h~ component gives 
the same result, and does not have to be considered separately. We can write this as a path 
integral expression: 



= (h+(y)) = J DiP+V?p-h + (y)e iSint[lpa] - iSint[lp e ] 

=Jv^+V^h+(y) [l + i |d 4 x(£+ -£f nt )-i |d 4 x(A + nt-A-„t) /dV(£ 1 + nt -£ i ; it ) + . 

= -ij d 4 xD+ + (y - x)A(x), (110) 

with if) a running over all fields. The equations of motion are then 

A(x) = Q. (Ill) 

The relevant interactions are given in Sec. I2.2I The one-loop equations of motion can be 
calculated order by order in perturbation theory, that is, ordered by the number of insertions 
coming from C- mt . The higher order expectation values can be evaluated by taking all possible 
Wick contractions. 



First order. At zeroth order there is no contribution because {h + ) = in the free theory. 
At first order there is a classical and quantum contribution. Starting with the classical tree- 
level contribution: 







-i J d 4 x(h + (y)h + (x))\+(x) = -i J d 4 xD+ + (y - x)X+(x) => A cl = X+(x) = 0. 



(112) 

The first order one-loop quantum contribution has only diagonal two-point insertions; the 
off-diagonal two-point interactions only enter at higher order. Then 

= -i J d 4 x{h + (y)h + (x)ifj+(x) 2 )^ml = -i J d 4 xD+ + {y - x)D+ + (0)^m 2 a 

A 1 = ±d (t> mlD++(0), (113) 

as there is only one possible Wick contraction. Here the subscript on A n denotes the nth 
order contribution. 



Second order. For convenience, in the rest of this section we will drop the ± superscript. 
It should be kept in mind that at all times one should sum over all possibilities, and the 
incoming h propagator is always D^ + (y — x). This summation is done explicitly in the main 
text. 

Consider the second order contribution, with one two-point insertion. The three-point 
vertex connecting to the incoming h field can either be in the first or second factor of C mt , 
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which cancels the 1/2! in front (symmetry under x f-> x'). Consider first two diagonal two- 
point insertions, and the three-point vertex diagonal as well: 



1 



% I d 4 xD h (y-x) 



d$m 2 a (x) / & A x'm 2 p(x')D a p{x - x')D a/3 (x' - x) 



(114) 



since there are two possible Wick contractions. The part between the square brackets is j42 iag . 
For the Minkowski contribution with w& A we get 

d 4 xd 4 x'2^ (x')(/ife) [(Or - H(T))h(x)] A (x)9(x)A Q (x')9(x')) (115) 



d 4 x[(d T - U(T))D h (y - x)] 



i2g l d x'm eo (x')D 00 (x - x')D ee {x - x) 



-i j d 4 xD h (y - x) 



i2g(-d T — %{t)) / d x'm go (x')D 00 (x - x')D ee (x' - x) 



as there is only one possible Wick contraction. To get the last expression we integrated by 
parts. The 6Ai d l h- vertex does not contribute (even at higher order and finite terms), as going 
through the same steps, we get a result of the schematic form d l A{t) = 0. 
Finally, for the diagram with one (m 2 ) * insertion we get 



~ j d 4 xd 4 x / ^m^(x)(m 2 ) 0i (x / )(^(x) 2 ^o(x , )^(^)) 

-i J d A xD h {y-x) -ifymj^x) j dV(m 2 ) 0i (x / ) J D At o(x - x^D^x' - x) 



(116) 



as there are two possible Wick contractions. Both Wick contractions give the same result. 

Third order. At third order the diagram with two mixed mass (m 2 )° l contributes to the 
divergent part in A. There is a symmetry under the interchange of the positions of the three 
vertices, (x,x',x"), and taking a definite order removes the factor 1/3! in (jllOp . The result is 



i J d 4 xdVdV'i<9 m 2 M (x)(m 2 )°XxO(mT^ 



-i Jd 4 xD h (y-x) 

^m 2 M (x) f d 4 x'd 4 x"m 2 0i (x')m 2 0l (x") ^ D w {x-x)D^(x'-x')D^(x"-x) 



(117) 

where the sum is over the four possible Wick contractions (|54p . The term between brackets 
is Af 338 . 
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